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L’étude de la mécanique des fluides remonte au moins à l’époque de la Grèce antique 
avec Archimède qui a découvert la notion de la poussée d’Archimède pour un fluide au repos. 
Aujourd’hui, la mécanique des fluides est l’un des domaines de la recherche les plus actifs avec 
de nombreux problèmes non résolus ou partiellement résolus, comme les problèmes de la 
pollution atmosphérique. La résolution des équations régissant ces problèmes complexes fait 
appel aux méthodes de résolution numériques. 

Comme tout problème de mécanique, la résolution d'un problème de mécanique des fluides passe 
par la définition du système matériel S, particules de fluide à l'intérieur d'une surface fermée 
limitant S. À ce système, on applique les principes et théorèmes généraux de mécanique et de 
thermodynamique à savoir : 

- principe de la conservation de la masse, 

- principe fondamental de la dynamique, 

- principe de la conservation de l'énergie. 

La mécanique des fluides concerne l’étude du comportement des fluides et des forces 
internes associées. 

Elle se divise en statique des fluides et dynamique des fluides : 

Statique des fluides : C’est l’étude des fluides au repos 
Dynamique des fluides : C’est l’étude des fluides en mouvement 

Un fluide peut être considéré comme étant formé d'un grand nombre de particules 
matérielles, très petites et libres de se déplacer les unes par rapport aux autres. Un fluide est donc 
un milieu matériel continu, déformable, sans rigidité et qui peut s'écouler. Parmi les fluides, on 
fait souvent la distinction entre liquides et gaz. 

Les liquides et gaz habituellement étudiés sont isotropes, mobiles et visqueux. La propriété 
physique qui permet de faire la différence entre les deux est la compressibilité. 

- l'isotropie assure que les propriétés sont identiques dans toutes les directions de l'espace. 

- la mobilité fait qu'ils n'ont pas de forme propre et qu'ils prennent la forme du récipient qui les 
contient. 

- la viscosité caractérise le fait que tout changement de forme d’un fluide réel s'accompagne 

d'une résistance (frottements). 
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I- Pression et Force 

Les particules qui forment un fluide ne sont pas immobiles les unes par rapport aux 
autres. Elles sont agitées de façon désordonnée ce qui provoque de nombreux chocs entre elles et 
avec les parois. Par ces chocs, le fluide applique une force sur les parois. Ces forces sont 
appelées forces de pression. 



S 


Considérons la figure ci-dessus représentant une enceinte contenant un fluide. Ce fluide 
exerce donc des forces sur chacune des parois. Ces forces sont dirigées vers l’extérieur de 
l’enceinte et sont perpendiculaires aux parois. Si on considère la face hachurée de surface S, le 
fluide lui applique une force F. On peut ainsi définir la pression P du fluide comme le rapport de 
cette force F et de la surface S : 


P = 


F_ 

S 


La pression représente donc la force qui s’exerce sur chaque unité de surface. 


II- Unités de pression 

Dans le système international, l’unité pour la force est le Newton, noté N, et l’unité de la 
surface est le mètre carré (m ). Par conséquent, l’unité dans le système international pour la 
pression sera le N/nr. 

Une unité a été inventée pour la pression : C’est le Pascal, noté Pa. On a par conséquent 
l’équivalence suivante : 1 Pa = 1 N/m 2 

Il existe de nombreuses autres unités de pression couramment employées dans l’industrie, 
bar : lbar=10 5 Pa 

Atmosphère (atm) : 1 atm=101325 Pa 

Millimètre de mercure (mmHg) ou torr : 1 mmHg=133,32 Pa 
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III- Notions de pressions absolue et relative 

Nous vivons dans un monde qui est baigné au sein d’un fluide : l’air. 

On désigne par pression atmosphérique la valeur de la pression de l’air ambiant. 

Cette valeur (que l’on mesure à l’aide d’un baromètre) fluctue en fonction des conditions 
météorologiques et de la zone géographique. Toutefois, la valeur de la pression atmosphérique 
oscille autour d’une valeur moyenne qu’on appelle pression atmosphérique normale qui vaut 
101325 Pa. 

Lorsque la pression d’un fluide est supérieure à la pression atmosphérique on dit que ce fluide 
est sous pression. 

Lorsque la pression du fluide est inférieure à la pression atmosphérique, on dit que le fluide est 
sous vide. 

Une pression nulle (P= 0 Pa) correspond à un vide parfait qui correspond en fait à une absence 
totale de particules (atomes ou molécules). 

On définit la pression relative, que l’on note F, par : P = P - P atm 

P est la pression absolue. P et P' sont toutes deux des pressions et ont par conséquent la même 
unité. 

IV- L'équation de l'hydrostatique 

IV-1. Formulation mathématique 

La pression varie avec la hauteur dans le liquide. Un point du fluide sera donc représenté 
par son altitude notée z. L’altitude est la coordonnée du point sur un axe vertical et dirigé vers le 
haut. On fixera de manière arbitraire l’altitude 0 (l’origine) sur cet axe. On veillera tout de même 
à choisir une origine pratique comme par exemple le fond d’un réservoir, le centre d’une pompe, 
etc. 

Pour exprimer à l’aide d’une relation mathématique l’évolution de la pression au sein d’un fluide 
au repos, il convient de respecter scrupuleusement les hypothèses suivantes : 

• Le fluide doit être au repos. 

• Le fluide doit être homogène. On ne peut écrire de relation qu’au sein d’un seul et même 
liquide. 

On considère au sein de ce fluide homogène et au repos, deux points distincts 1 et 2, d’altitudes 
respectives zi et zi, alors on peut écrire la relation suivante entre les pressions P i et /L : 

^*1 P fluide •<?•-! — ^2 + P fluide S -Z 2 
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1 * 


1 * 


L x 



Où : g est l’accélération de la pesanteur. 

Cette relation est appelée équation de l'hydrostatique. On peut aussi l’écrire de la façon 
suivante : 

P + p.g. z = constante 


IV-2. Exemple 

La pression de l’eau dans l’océan à 15 m de profondeur est de 4 bar. On souhaite 
connaître la pression qui règne à 200 m de profondeur. La densité de l’eau de mer est de 1,02. 



A * 


B . 

On choisit l’origine des altitudes à la surface de l’océan. Le point A correspondant à une 
profondeur de 15 m aura donc une altitude négative Za=-15 m. De même, zb=- 200 m. On aura de 
plus : 

Pmer=d mer .p eau =1,02*1000 = 1020 kg/nP 
P A = A bar = 400000 P a 

On peut raisonnablement considérer l’océan comme un fluide homogène et au repos (on 
supposera qu’il n’existe pas de courants marins). Par conséquent, on peut appliquer l’équation de 
l’hydrostatique : 

P A+ Pmer-8-Z A =P B +P mer -g-Z B 
D’où : 

p b — p a +Pm er -g{z A -zj = 400000 + 1020 *9.81 *(-15- (-200)) 

P R = 225 1 147 P =22.5 bar 

jts a 
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La pression de l’océan à 200 m de profondeur est donc de 22.5 bar. 


IV-3. Remarque : cas des gaz 

Les masses volumiques des gaz sont très faibles en comparaison de celles des liquides. 
Par conséquent, dans un gaz, la pression variera très peu avec l’altitude. On pourra donc 
considérer en génie des procédés que dans un gaz la pression est partout la même. 

IV-4. Interface entre deux fluides 

Il est très courant de rencontrer deux fluides en contact comme un gaz au-dessus d’un 
liquide ou un liquide léger surnageant sur un liquide plus lourd avec lequel il n’est pas miscible. 
Cette surface de contact est, pour un liquide au repos, plane et horizontale. On l’appelle aussi 
interface. 

On aura alors la propriété suivante : la pression est identique de part et d’autre de l’interface. En 
effet, lorsque l’on traverse l’interface, on n’observe pas de discontinuité brutale de la pression. 
Ainsi, par exemple, la pression d’un liquide à sa surface ouverte à l’atmosphère est la pression 
atmosphérique. 



Si on considère la figure suivante, on peut voir l’interface entre le fluide 1 et le fluide 2. Le point 
A est situé à l’interface des deux fluides. Par conséquent, le point A appartient à la fois au fluide 
1 et au fluide 2. La pression P A est donc la même que l’on considère que le point A fasse partie 
du fluide 1 ou du fluide 2. On peut ainsi appliquer l’équation de l’hydrostatique au sein du fluide 
1 entre les points A et B, et aussi sur le fluide 2 entre les points A et C : 

P A — P B P\ -8-{Zb ~ Z A ) 

P A = P C +p2-g{z c -Z A ) 

Par conséquent, on aura la relation suivante entre P b et Pc : 
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P B = P c + Pi -S-{z c ~ Z a)~ P\ ■ g.{z B - Z A ) 

On notera enfin la remarque suivante qui est très importante : on ne peut pas écrire directement 
l’équation de l’hydrostatique entre les points B et C, car ils appartiennent à des fluides 
différents. 


V- Corps immergés : Principe d'Archimède 

Analysons ce qui se passe du point de vue de l’équilibre mécanique sur un flotteur, c’est à 
dire réfléchissons aux forces auxquelles il est soumis. 


Poussée 

d’Archimède 


Flotteur 



Poids 


Partie 
immergée 
{volume V in 


Liquide 


Tout d’abord, le flotteur est soumis à son propre poids, dirigé vers le bas et d’intensité m ' g, où 
m est la masse du flotteur et g l’accélération de la pesanteur. 

Mais le poids n’est pas la seule force, sinon le flotteur coulerait. Il est donc soumis à une force 
dirigée en sens inverse du poids et égale en intensité pour assurer l’équilibre des forces. 

Cette force est appelée Poussée d’Archimède. Le principe d’Archimède permet de relier 
l’intensité de cette poussée au volume de flotteur immergé et à la masse volumique du liquide. 
L’énoncé traditionnel du principe d’Archimède est le suivant : 

Tout corps plongé dans un liquide subit une poussée verticale dirigée vers le haut dont 
l’intensité est égale au poids de liquide déplacé. 

Il est toutefois possible d’en donner une présentation plus simple. Ainsi, si l’on appelle V imm , le 
volume de la partie du flotteur qui est immergée dans le liquide, alors l’intensité A de la poussée 
d’Archimède s’exprimera de la façon suivante : 

P P liquide ^ imm 
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Remarque : 

Pour que le flotteur soit en équilibre, c’est à dire pour qu’il flotte, il faut que la poussée 
d’Archimède soit égale en intensité au poids du flotteur. Or, le poids du flotteur vaut : 

m flotteur S- 

On aura donc : m flotteur .g = p liqmde .g.V imm 

En remplaçant m flotteur par p flotteur .V flone!tr , on obtient : 

P flotteur ^ flotteur P liquide ‘ Q ^ imm 

Le volume immergé V imm est toujours plus petit que le volume total du flotteur V flotteur . Par 
conséquent, pour que la précédente relation soit satisfaite et que le flotteur flotte, il faut que la 
masse volumique du liquide soit plus faible que celle du liquide, p flotteur < p liquide 

On pourra encore écrire cette condition avec les densités : d flotteur < d Uquide 


VI- Action d’un fluide sur une paroi 

VI-1. Force exercée sur un élément de surface 

Soit un liquide contenu dans un récipient de forme quelconque. Soit P 0 la pression du gaz 
environnant le récipient. 

En un point M situé à une distance h de la surface libre du liquide, la pression est : 

P(M) = P ü +pgh 

L’élément de surface dS entourant lepoint M subit : 

Une force (lf\ normale à dS et orientée vers 
l’extérieur de la part du liquide : 

df\ = PdS =[P 0 +pg h] dS 

Une force df 2 normale à dS et orientée vers 
l’intérieur du fluide de la part du gaz environnant 

df 2 = P Q dS 

La résultante s’exerçant sur dS est donc [force relative] : 

df = p g hdS soit df = p g h dS n 
Où : n est un vecteur unitaire normal à dS et orienté vers l’extérieur. 
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VI-2. Cas d’une paroi plane 

Soit l la largeur de la paroi perpendiculaire au plan 
de la figure et L sa longueur dans le plan de la figure. 

On a : 


dS = dl dL 

dz 

dL 


sina ■ 


dL 


dz 

siner 


Or: 



z = L sin a 


dS = dl 


dz 

sin a 


df - p g zdl — 


dz 


Donc 


F = 

P g f' 


sin a •’° 

F = 

PS l: 


sin a 


sin a 


l pLs'ma 


[ dl f 

Jo Je 



dz 


0 Jo 
2 2 


z dz 


3111 , .1 13111 Wj 

= p g l L~ 

2 2 

Détermination du point d’application 

Toutes les forces élémentaires sont orientées suivant n 


— , T n sina 

F = p g IL- — — 


la résultante l’est aussi : 


Pour déterminer le point d’application de F , appelé centre de poussée C, on écrit le moment en O 
ÔC aF = ^ÔP Adf 

Le moment de la résultante = résultantes des moments en un point (O par exemple). 

ÔC= OC u 
OC • F = J OP ■ df 


OC-pglL- 


OC p g l L 2 


OC = 


sina 


-- . 

P g z 


J sin a 


sin a 

, Pg 

L . 2 

rLsina 

z 2 dz 

sin" a 

Jo 


rL sin a 



i t 

z 1 dz 



sin 2 a 
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=> oc 


1 Ü sin 3 a 
sin 2 <2 3 


sin a x2 

L 

2 

2 L 

=> OC = -y- Indépendant de a 

VI-3. Cas d’une paroi Cylindrique 



On considère le schéma ci-dessus. 

Où : Oe- est l’axe perpendiculaire au plan de la figure et 1 la largeur suivant cet axe. 

K 

L’angle 9 varie de 0 à 

Le vecteur normal à la paroi est : n — e r 

Toutes les forces relatives élémentaires passent par O, d’où la résultante est appliquée 
en O = C centre de poussée. 


clf = p g h dS n 
Avec : dS = R dû dz 

En prpjetant cette relation suivant la verticale et l’horizontale, on aura : 
n = e r = cosû i + sin# j 

h 

Or : h = R cos 0 car cos 0 = — 

R 


— > O centre de poussée ; en effet : 

ÔCAF=^ÔMAdf = 0 


Or 


OM // df VM g paroi => O = C 
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— » Composante horizontale F y 

df = p g h dS e r • e = p g h dS sin 9 
F y = J df y = p g J R cos 9 R dû dz sin 9 

ci 

F y = p g R 2 ly 2 cos 0 sin 9 dd = p g R~ / — 

— > Composante verticale F x 

F x = J df x = avec df x = p g h dS cos 0 


F x = p g J R cos 6 R dd dz cos 9 = p g R 1 /j* 2 cos 2 9 dd 


Or: 

cos 2 # + sin 2 0=1 
cos 2 0 - sin 2 # =cos2 9 
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I- Description d’un fluide en mouvement 


1-1. Représentation eulérienne et représentation lagrangienne 
1-1.1 Représentation lagrangienne 

Décrire le mouvement d’un fluide fait appel à des notions différentes de celles 
développées en Mécanique du point ou du solide. Le mouvement d’un fluide est un écoulement 
où il y a déformation continue du fluide. 

On peut, de manière analogue à ce que l’on fait en Mécanique du solide, isoler (par la pensée ou 
en trouvant un moyen de visualisation, coloration par exemple) une partie restreinte du fluide 
appelée particule et la " suivre " au cours du temps c’est à dire connaître à chaque instant sa 
position. 

Cette position sera connue, par exemple, par ses coordonnées cartésiennes : 

x p (t,x 0 p ,y 0 p ,z 0p ), y p (t,x 0 p ,y 0 p ,z 0p ) et z p (t,x 0 p ,y 0 p ,z 0p ) 


Où : x 0 ,y Q ,z 0 représentent les coordonnées de la particule choisie à l’instant t 0 , la vitesse de 


la particule aura pour composantes : v x 


d - x P 

"aT’ Vy 


dt 


■et v_ 


dt 


Au cours du temps, la particule sera en différents points M, l’ensemble des points M constitue la 
trajectoire de la particule. 

Cette façon de faire est appelée méthode de Lagrange, les variables introduites sont appelées 

variables de Lagrange. 

La méthode de Lagrange s’avère dans la plupart des cas délicate, car il n’est pas facile de suivre 
les particules : elle est peu employée. 

1-1.2 Représentation eulérienne 

La méthode d’Euler consiste à connaître la vitesse des particules au cours du temps un 
endroit donné déterminé par ses coordonnées, par exemple cartésiennes x, y, z. Elle est plus 
employée que la méthode de Lagrange, la connaissance du champ des vitesses étant suffisante 
pour la description du fluide en mouvement. 

Les composantes du vecteur vitesse v sont des fonctions des variables (x, y, z), ainsi 
v = v x i + v y j+v.k 

: v x =v x {x,y,z,t ) , v y = v y (x,y,z,t) et v, =v z {x,y,z,t) . 

Remarque : 

Dans la méthode d’Euler, l’accélération d’une particule peut être due, bien sur, au caractère 
instationnaire de l’écoulement, mais aussi à sa non-uniformité. Ainsi, chacun a pu constater, dans 
l’écoulement permanent d’une rivière, l’accélération des particules lors du franchissement d’un 
rétrécissement. 


20 


II- Définitions 


II- 1. Ecoulement permanent 

L’écoulement du fluide est permanent ou stationnaire si ses composantes de vitesse 
sont indépendantes de la variable temps t ; il est dit non-permanent ou instationnaire si cette 
condition n’est pas réalisée. 

II-2. Ecoulement uniforme 


L’écoulement du fluide est uniforme si ses composantes de vitesse sont indépendantes 
des coordonnées d’espace ; il est non-uniforme si cette condition n’est pas remplie. 

II-3. Ecoulement plan 


Un écoulement plan est un écoulement dont le champ de vitesse est, à tout instant, 
parallèle à un même plan, et qui ne varie pas perpendiculairement au plan. 

v, {x l ,x 2 ,t) 

On a donc : v(x,t) = \v 2 (x, , x 2 , î ) 


II-4. Ecoulement irrotationnel 

Un écoulement est dit irrotationnel si rot(v ) = 0 . 


Dans ce cas, 3(p appelée potentiel des vitesses telle que v = - grad{(p ) 

Remarque : 


Le vecteur cb = ^rot(v) est appelé vecteur tourbillon. 


II-5. Ligne de courant à un instant t Q fixe 

On appelle ligne de courant une courbe dont la direction tangente en chacun de ses 
points est la direction du vecteur vitesse. L’équation d’une ligne de courant se calcule par 
intégration des équations obtenues à l’aide de : 



v. 

clx 

VMM = 

V 2 A 

dy = 


V 3 

dz 


II- 6. Tube de courant 
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Un tube de courant est un ensemble de lignes de courant s’appuyant sur un contour 

fermé. 


II-7. Ligne d’émission 


instant donné par des particules passées antérieurement en un même point. 
Trajectoire, ligne de courant et ligne d’émission sont confondues pour un écoulement permanent. 

II-8. Trajectoires 

Ensemble des positions x occupées par une particule fluide donnée. 
dix 

Elle est donc solution de : — = v(x,t) <=> dx j = v j (x,t)dt 

dt 

Ce qui conduit à (sivj ^ 0,v 2 ^ 0 etv 3 ^ 0) : 


On a trois équations du premier ordre, donc 3 constantes d'intégration. On obtient ainsi une 
famille de courbes à 3 paramètres. 

Pour observer, au sens propre, des trajectoires, on peut mettre en suspension dans le milieu 
quelques particules et faire une photographie avec un temps de pose très long. 

Remarque : 

A priori, les trajectoires et les lignes de courant sont des entités différentes. Toutefois ces deux 
concepts sont identiques dans le cas d'écoulements stationnaires. 

III- Equation de continuité 

III-l. Dérivation suivant la méthode d’Euler 

Considérons la fonction scalaire Q(x, y, z, t) rendant compte d’une grandeur physique 
caractéristique du fluide au point de coordonnées x, y, z et au temps t. 

La particule fluide au temps t + dt sera au point de coordonnées x + v x dt, y + v dt,z + v.dt . 
La variation de la fonction Q sera donc égale à : 


On appelle ligne d’émission une courbe constituée par l’ensemble des points atteints à un 


dx j _ dx 2 _ dx 3 


v x {x,t) v 2 (x,t) v 3 (x,t) 



dx dy dz 


v.dt + —v.dt + — dt 
y a, z 


9g 


dt 



Cette dérivée apparaît comme la somme de deux termes : 
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le premier, qualifié de convectif ou advectif, est du à la non-uniformité de l’écoulement, 
le second, qualifié de temporel, est du au caractère instationnaire de l’écoulement. 


III-2. Théorème de la divergence 

Nous nous contenterons de donner, sans démonstration, un énoncé de ce théorème appelé 
aussi théorème de Green - Ostrogradski : 

Le flux d’un champ tensoriel A au travers de la surface dD enveloppant le domaine D est 
égal à l’intégrale de la divergence du champ tensoriel sur le domaine : 

J" An ds = ^divëkdv 
as s 


Remarque : 

Dans le cas où A représente un champ vectoriel constant, on obtient J nds 


dD 


Dans le cas où A représente un champ scalaire f, on obtient J f nds = J gmd(f)dv 

a d d 

III-3. Théorème de l’intégrale nulle 

L’énoncé de ce théorème est le suivant : 

Considérons un champ vectoriel volumique a défini et continu sur un domaine D. Si 
quelque soit le sous domaine D* inclus dans D, l’intégrale du champ tensoriel sur le 
domaine D* est nulle, alors le champ tensoriel est identiquement nul. 

Jllüfv = 0 VD' a=§ 

III-4. Equation de continuité 

Le principe de conservation de la masse postule qu’il n’y a ni apparition ni disparition de 
matière. En conséquence la variation de la masse au cour du temps est nulle : 


dM 

di 


= 0 


La masse peut se calculer à partir de la masse volumique : 

d 




ûdv = 0 


23 


Avec la notion de dérivée particulaire d’une intégrale de volume, on obtient : 


JVüfv = dv + üfe = 0 


On peut encore utiliser le théorème de la divergence : 

j a a 

— | | dfv + ds = | ûfv + |üfîV'(,i7V) ifv 


= 0 


Ce qui nous donne une forme locale de l’équation de continuité avec le théorème de l’intégrale 
nulle : 

— —+div(, 0 v)=O 
dt 


De plus, nous avons les relations : 


div{ ^v) = ,û<div(v) + v . grad ,& 
d ,0 3,tf _ 


= — —+v.grad ,o 


dt dt 


On obtient ainsi une autre forme locale de l’équation de continuité : 
d ,0 


dt 


■+ ,üdiv(y ) = 0 


Remarque : 

Sauf précision contraire, nous appliquerons l’équation de conservation de masse en absence de 

source ou de puits, soit — + pdiv(pv) = 0 

dt 

Deux cas particuliers sont alors à considérer. 

Le cas 1 d’un fluide incompressible (p = cte) => div(v) = 0 pour un écoulement stationnaire 
ou instationnaire. 

Cet écoulement est dit isovolume. 

Le cas 2 d’un écoulement stationnaire — = 0 : 


dt 


div(pv) = 0 = p div(v ) + v grad (p) 


En dehors de la possibilité cas 1, il existe la possibilité d’écoulements isovolumes tels que 
v grad(p) = 0 où les variations de masse volumique sont orthogonales, en tout point, au vecteur 
vitesse. 

Ce cas correspond à des écoulements stratifiés par salinité ou température (courants marins), par 
température et humidité (atmosphère). 
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III-5. Vitesse et accélération 


III-5.1 Vitesse 


Pour une particule fluide donnée, identique pour M fixé, la vitesse est donnée par 

v[m ,t) = 


Soit en notation tensorielle : 


v 


dx 
dt 


M 


M, 


III-5.2 Accélération 


Pour une particule fluide donnée, identique pour M fixé, l’accélération est donnée par 




dv 

dt 


Soit en notation tensorielle : 


M 

d\’: 


r t) 


III-6. Conservation du débit 


111-6.1 Définitions 

i. Le débit 

Le débit est le quotient de la quantité de fluide qui traverse une section droite de la conduite par 
la durée de cet écoulement. 


ii. Débit-masse 

Si dm est la masse de fluide qui a traversé une section droite de la conduite pendant le temps dt, 


par définition le débit-masse est : q m 
Son unité est le kg-s' 1 


dm 

dt 


iii. Débit-volume 

Si dV est le volume de fluide qui a traversé une section droite de la conduite pendant le temps dt, 


par définition le débit-volume est : q v = 


dV 

dt 
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3 1 

Son unité est le m -s' . 


vi. Relation entre q m et qv 

dm 

La masse volumique est donnée par la relation : p = d'où: q m =pq v 

dV 

Remarque : 

Les liquides sont incompressibles et peu dilatables (masse volumique constante) ; on parle 
alors d'écoulements isovolumes. 

Pour les gaz, la masse volumique dépend de la température et de la pression. Pour des vitesses 
faibles (variation de pression limitée) et pour des températures constantes on retrouve le cas d'un 
écoulement isovolume. 



III-6.2 Conservation du débit 

Considérons un tube de courant entre deux sections Si et Si. Pendant l'intervalle de temps 
dt, infiniment petit, la masse dmi de fluide ayant traversé la section Si est la même que la masse 
dm 2 ayant traversé la section S 2 . 

En régime stationnaire, le débit-masse est le même à travers toutes les sections droites d'un 

même tube de courant. q ml = q m2 

Dans le cas d'un écoulement isovolume (p = Cte) : 

En régime stationnaire, le débit-volume est le même à travers toutes les sections droites d'un 
même tube de courant. q vx = q v2 

III-6.3 Expression du débit en fonction de la vitesse v 

Le débit-volume est aussi la quantité de liquide occupant un volume cylindrique de base 
S et de longueur égale à v, correspondant à la longueur du trajet effectué pendant l'unité de 
temps, par une particule de fluide traversant S. 
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Il en résulte la relation importante : q v =V S 


III-6.4 Vitesse moyenne 




If tfmcy 


En général, la vitesse v n'est pas constante sur la section S d'un tube de courant ; on dit qu'il 
existe un profil de vitesse (à cause des forces de frottement). Le débit-masse ou le débit-volume 
s'obtient en intégrant la relation précédente : 

Dans une section droite S de la canalisation, on appelle vitesse moyenne V la vitesse telle 

q ue : v ,noy = Y 

La vitesse moyenne V moY apparaît comme la vitesse uniforme à travers la section S qui assurerait 
le même débit que la répartition réelle des vitesses. 

Si l'écoulement est isovolume, cette vitesse moyenne est inversement proportionnelle à l'aire de 

la section droite. q v =V hnoy S l =V 2moy S 2 =cte 

C'est l'équation de continuité qu’on peut aussi écrire sous la forme : 

^2 

La vitesse moyenne est d'autant plus grande que la section est faible. 


III-7. Théorème de BERNOULLI 
III-7.1 Le phénomène 
Observations 

- Une balle de ping-pong peut rester en suspension dans un jet d'air incliné. 

- Une feuille de papier est aspirée lorsqu'on souffle dessus. 

Conclusion : La pression d'un fluide diminue lorsque sa vitesse augmente. 
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III-7.2 Théorème de Bernoulli pour un écoulement permanent d’un fluide parfait 
incompressible 

Un fluide parfait est un fluide dont l'écoulement se fait sans frottement. 

On considère un écoulement permanent isovolume d’un fluide parfait, entre les sections Si et S 2 , 
entre lesquelles il n’y a aucune machine hydraulique, (pas de pompe, ni de turbine). 

Soit m la masse et V le volume du fluide qui passe à travers la section Si entre les instants t et 
t+dt. Pendant ce temps la même masse et le même volume de fluide passe à travers la section S 2 . 
Tout se passe comme si ce fluide était passé de la position (1) à la position (2). 

En appliquant le théorème de l’énergie cinétique à ce fluide entre les instants t et t+dt (la 
variation d’énergie cinétique est égale à la somme des travaux des forces extérieures : poids et 
forces pressantes), on obtient : 

y 2 

P — + pgz + P = Cte 
V 2 

p est la pression statique , pgz est la pression de pesanteur , p— est la pression cinétique . 

Tous les tenues s’expriment en pascal. 

En divisant tous les termes de la relation précédente par le produit g, on écrit tous les termes dans 
la dimension d'une hauteur (pressions exprimées en mètres de colonne de fluide). 

P V 2 

H est la Hauteur totale, — est la Hauteur de Pression, z est la cote, — est la Hauteur 

Pg 2g 

P 

cinétique , z est la Hauteur piézométrique . 

Pg 



III-7.3 Démonstration : 

On considère le schéma suivant : 
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k- V-2&1 — *2 -* 



L'équation de Bernoulli pour les fluides incompressibles peut être démontrée, soit, par 
application du principe de conservation de l'énergie le long d'une ligne de courant, en négligeant 
les effets thermiques, de viscosité, de compressibilité, soit, à l’aide de l’équation d’Euler. 

i. A l’aide de la conservation d’énergie 

Soit le système fermé contenu à l'instant t entre x\ et %2 et à t + At entre x\ + v\ .At et X2 + 

V 2 -A t. 

Le fluide est incompressible, la masse Am contenue entre xi et xi doit être identique à la masse 
contenue entre x\ + v\.At et X2 + V2.A t. 

Ce que l'on peut ramener ici à la conservation du débit massique : 

At.v1.A1. p = At.v2.A2-p. 

Toutes les forces qui s'exercent (forces pressantes et poids) sont conservatives (il n'y a pas d'effet 
visqueux). On peut donc appliquer le théorème de conservation de l’énergie mécanique au 
système: 


AE PP + A E k =W 
Où 

— Ani( t'9 — t'f ) / - e st la variation d’énergie cinétique du système. 

A E pp = Amgh2 - Amghi est la variation d’énergie potentielle de pesanteur du système. 
W = p\A\.{viAt) - P2A2A V2A1) est le travail des forces de pressions. 

Soit: 


-Amvn 
2 2 


mvl + A mgh 2 - Amghx = p l Ai.(v i At) - p 2 A 2 .(v 2 At) 
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D'où, en divisant par Am: 


»1 , . ,P'- V î, nh , 

— + H — TT + + 

1 pl p 


Et donc: 


V V 

+ g h + L = C 

2 p 

On peut remarquer que la démonstration est faite dans le contexte particulier d'un écoulement 
obéissant à la géométrie de la figure. Cependant, pour un écoulement quelconque en régime 
permanent, on pourra toujours définir au voisinage d'une ligne de courant une section sur 
laquelle la vitesse est homogène, et raisonner comme précédemment. 


ii. A l’aide de l’équation d’Euler 

Quelque soit la particule fluide M, l’équation d’Euler peut s’écrire sous la forme : 
pf{M / R) = -gracl(P + pgz) 

Où : z est la côte de la particule M et R un référentiel galiléen. 

D’autre part, on sait que : 


y(M / R)= grad 


r v 2 ^ 

v^y 


+ rot 


(v)aV 


Dans le cas d’un écoulement permanent, p = Cte, ce qui conduit à : 


grad 


( y 2\ 

v ^ y 


+ p rot {y ) a V = —grad (P + pgz) 


Comme la fonction grad est linaire, alors on aura : 


grad 


f V 2 

p— + P + pgz 


+ prot 


(y) a y = o 


Considérant un déplacement dM infinitésimal sur une ligne de courant (ou sur les trajectoires, 
puisque l’écoulement est permanent, les trajectoires et lignes de courant sont confondus), on a : 

dM = Vdt 

Comme le vecteur vitesse V est parallèle à l’élément de déplacement dM et que V est 


perpendiculaire à rot {y) aV , alors la projection de grad 
suivant le vecteur dM conduit à: 


f y 2 

p—+P + pgz 


+ prot 


(y) a y = o 
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grad 


f y 2 ) _ 

p — + P + pgz dM =0 

V ^ y 


On sait que quelque soit la fonction scalaire f, on peut écrire : 


grad(f)dM -df 


Alors, on peut déduire que : 


f y 2 


p— + P + pgz 


= 0 Sur la ligne de courant considérée 

v z ; 

Et par conséquent : 

' V 2 ) 

p P + pgz = Cte Sur la ligne de courant considérée 

V ^ 

La constante varie d’une ligne de courant à une autre. 


III-7.4 Equation transversale 

Soit N le vecteur normal à la trajectoire (ou à la ligne de courant) orienté vers l’intérieur 
de la cavité. 



Dans le repère de Frénet (f , n), on a : 

/ x dV - V 2 - 

y(M/R) = T + — N 

dt R 

Où : R est le rayon de courbure des lignes de courant au point M considéré. 


La projection de p f = -grad (p) - pge. suivant la normale conduit à : 

V 2 dP / X 

p— = - — -pgcos{(p) 

R an 


Où : (p est l’angle entre les vecteurs unitaire N et <?_ . 
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Remarque : 

Dans un écoulement où on peut négliger les forces de gravite, la pression reste constante 
lorsqu’on se déplace perpendiculairement aux trajectoires dans les zones ou celles-ci sont des 
droites. C’est-à-dire lorsque R tend vers l’infini. 

En général, si les trajectoires sont des droites alors on aura : 


III-8. Applications 

III-8.1 Phénomène de Venturi 

Un conduit de section principale S A subit un étranglement en B où sa section est S B . La 
vitesse d’un fluide augmente dans l’étranglement, donc sa pression y diminue : 
si v B > v A donc p B < p A 
Le théorème de Bernoulli s'écrit ici : 

Pà +^pvi = P B PVb = Pc +^pv c 

D'après l'équation de continuité, V B S B =V A S A = q v et V B > V A donc P A > P B 

1,1 1 v 2 , 2 

pA-p B = 2p q = k q 

B A 

La différence de pression aux bornes aux extrémités du tube de Venturi est proportionnelle au 
carré du débit ; application à la mesure des débits (organes déprimogènes). 

On peut citer aussi la trompe à eau, le pulvérisateur... 



III.8.2 Écoulement d'un liquide contenu dans un réservoir - Théorème de Torricelli 

Considérons un réservoir muni d'un petit orifice à sa base, de section s et une ligne de 
courant partant de la surface au point (1) et arrivant à l'orifice au point (2). En appliquant le 
théorème de Bernoulli entre les points (1) et (2), 
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p^-ppgz, + p, = p^-+ pgz 2 + p 2 

Or : pi = p 2 = pression atmosphérique et vi«v 2 d'où : V 2 = gz 

La vitesse d'écoulement est la même que la vitesse de chute libre entre la surface libre et l'orifice, 
quelle que soit la masse volumique du liquide. 

Application : vase de Mariotte à débit constant. 


III-8.3 Vidange d’un réservoir 

On considère le schéma suivant : 


La section S = S(z) est variable. D’après, d’une part, la conservation du débit S.V = s.v et d’autre 
part, Torricelli v = ^2gz , on aura : 

SV = Sy[2gz 

dz 

Or : l’expression de la vitesse à la surface libre s’écrit : V = — 

dt a 

Donc, il vient : 


f dz} 

V dt j 


ZI 




Cette équation peut s’écrire sous la forme : 

Î + W )^= 0 


vi=0 


3 Vï 


jet 

parabolique 


La section S (z) peut s’écrire sous la forme générale : Bz n 
En particulier : 

le cas d’un réservoir parallélépipédique : S(z) = L*z 
le cas d’un réservoir conique : S(z) = 7Ttg 2 (a)z 2 


Où : a est l’angle d’ouverture du cône 


L’équation précédente peut s’écrire alors sous la forme : 

— + — 72g! = 0 

dt Bz n 

L’intégration de cette équation entre deux altitudes z, et z 2 correspondant aux instants /, et t 2 
conduit à : 


h ^2 


Remarque : 
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Si H est la hauteur initiale de la surface libre, alors la durée nécessaire pour vider le réservoir 



III-9. Théorème global d’Euler 

On considère un tube de courant au sein d’un fluide parfait incompressible de masse 
volumique p en écoulement permanant. 

La portion étudiée est délimitée par Si, S2 et Su son volume est V. 


Si 



Le volume de la tranche hachurée est 
dV = S dx = S dx = S n dx 

Or : dx = v dt => dV = S ( v - n) dt = Q dt 

Où : Q est le débit volumique 

Sa masse est donc : dm = p dx = p Q dt 

Sa quantité de mouvement : dp = ( dm)v = p Qv dt 

Le taux de variation de la quantité de mouvement est : — — = p Qv 

dt 

En utilisant la relation fondamentale de la dynamique (R.F.D) à la tranche hachurée, on aura : 

{dm) y{%)= ext appliquées %R 

(dm) ^V(%)= Yje, 
dt 
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Où : R est la résultante de tous les efforts appliqués 

Or: Q, = Q 2 => /?ô(v 2 -v 1 ) = R 

La résultante R est la somme des forces suivantes : 
+ forces de pression 

+ forces de pesanteur 
+ forces exércées par des obstacles 
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I- Introduction 

Le formalisme de Lagrange permet d’étudier une vaste gamme de problèmes en 
mécanique. Il est équivalent au formalisme de Newton mais, il a sur ce dernier un certain nombre 
d’avantages. D’abord, il est fondé sur un principe théorique fondamental et facile à retenir. Il 
utilise des quantités scalaires plutôt que vectorielles et sa forme est indépendante des 
coordonnées utilisées. 

On considère un système de N particules, numérotées (i = 1,2,..., N), de masse m, soumis à un 
ensemble de forces. Nous distinguons deux types de forces : externes et internes. 

> Les forces externes, agissant sur le système, sont dues à des sources externes au système. 

> Les forces internes, agissant sur la particule i, sont dues à toutes les autres particules du 

système. 

Ainsi l’équation du mouvement 

forme : — (m ; V. ) = F iext + ^ F jt 

Où, F iext représente une force extérieure agissant sur la particule i ; 

Et F jt représente une force interne agissant sur la particule i, due à la particule j. 

Evidement, nous supposons que F jt obéit à la loi de l’action et la réaction. C’est à dire que 
lorsque la particule j exerce une force F ji sur la particule i alors la particule i exerce une force 
F; j sur la particule j . tel que: F ij + F Ji = 0 

Dans un référentiel galiléen donné, il faut, pour décrire le mouvement du système, connaître la 
dépendance temporelle des 3N coordonnées r ; O,- , y t ,z t ) . 

Les lois de Newton permettent de décrire le mouvement de chaque particule localement et de 
proche en proche tout au cours de temps. Par contre, le formalisme de Lagrange permet de 
décrire le mouvement du système globalement. 

II- Coordonnées généralisées 

Afin de déterminer la position d’un système de N points matériels dans l’espace il faut 
donner N rayons vecteurs, c’est-à-dire 3N coordonnées. 

En général, le nombre de grandeurs indépendantes qu’il faut se donner pour déterminer de façon 
univoque la position d’un système est appelé nombre de degré de liberté du système. Dans 
d’autre cas, le nombre de degré de liberté est égal à 3N. Ces grandeurs ne sont pas forcement les 


pour la particule i aura la 

j * i 
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coordonnées intérieures du point, et selon les conditions et la nature du problème, le choix d’un 
autre système de coordonnées sera plus utile. 

Soit S grandeurs quelconques, qi, q2,....q s , qui caractérisent complètement la position du système 
(à S degré de liberté) sont appelées coordonnées généralisées et les dérivées ( q t i = l,s) ses 
vitesses généralisées. 

> Les coordonnées d’un point matériel sont utilisées pour repérer un point par rapport à un 
référentiel de référence. 

> Les coordonnées généralisées sont utilisées pour repérer un système dynamique par 
rapport à un référentiel de référence. 

Exemple de coordonnées généralisées : 


Remarque : 

Un ensemble de coordonnées généralisées est dit complet si les valeurs des coordonnées, 
correspondant à une configuration du système compatible avec les liaisons, sont suffisantes pour 
repérer la position de chacune des parties du système. 

III- Contraintes géométriques 

III- J. Définition 

Une contrainte géométrique est une condition qui entraîne une réduction du nombre de 
degrés de liberté. 

Exemple : 

> Roulement sans glissement d’un disque sur un plan. 

> Une tige assujettie à rester dans un plan pendant son mouvement. 



e 


y 


*- 
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III-2. Liaison holonome 

Une liaison est holonome si elle s’exprime par une relation entre les variables de position 
et le temps sans faire intervenir les dérivées par rapport au temps des variables de position. Une 
liaison holonome est donc de la forme : 

f(q l ,q 2 , -,q N ,t) = 0 


Exemple : 

On considère le mouvement d’un disque sur 
un axe horizontal représenté par le schéma 
suivant : 

On a : x = R 6 c’est une liaison holonome 
c’est-à-dire dx = Rd6 



III-3. Contrainte non holonome 

Une liaison est non holonome si elle s’exprime par une relation entre les variables de position, le 
temps et les dérivées premières, par rapport au temps, des variables de position. Une liaison non 
holonome est donc de la forme : 


Remarque : 

Une liaison est non holonome si elle n’est pas intégrable. 

IV- Centre de gravité 

IV -1. Cas d’un ensemble discret 

Le centre de gravité d’un ensemble discret de N particules, numérotées (1,2,..., AO et de masse 
est G défini par : 


OG = 


m j OG { + m 2 OG 2 H 1- m N OG N 

m l +m 2 -\ m N 


IV-2. Cas d’un ensemble continu 
Le centre de gravité G d’un ensemble continu est défini par : 
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OG = 



P G système 


V- Le paramétrage 

On commence par préciser les paramètres de rotation et ensuite ceux de la translation 

V-l. Paramètres de rotation 

V-l.l Rotation propre (cp) 

Définition : 

C’est la rotation du système autour de l’axe de révolution. 

Remarque : 

Lorsque le système est rectiligne (cas d’une tige) il n’y a pas de rotation propre. 

V-l. 2 La nutation (0) 

Définition : 

C’est la rotation du système soit autour de ü , soit autour de v . La nutation définit l’angle entre 
l’axe fixe z 0 et l’axe de révolution Z du système. 

Remarque : 


> Si la nutation est autour de u alors : 




> V 


z 0 = cos 6 z - sin 6 w (2) 


> Si la nutation est autour de v alors : 



► u 


w 
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Remarques : 


> Si le mouvement est plan alors la rotation (nutation) est égale à une constante 

> Si l’axe de révolution du système reste parallèle avec le plan fixe alors la 


nutation est égale à une constante 


V- 1.3 la Précession (\|/) 


C’est la rotation du système autour de l’axe fixe z 0 


Définition : 


v 



* *0 


u 


V-2. Paramètres de la translation 

On définit ces paramètres en exprimant le vecteur OG dans le référentiel 
fixe 

Remarque : 

Si le système admet un point fixe, on n’aura pas besoin de définir ces paramètres. 

V-3. Schéma du paramétrage 


Translation OG 



*- 



z V 


Z u 



Le vecteur rotation du système est le suivant : 



0 u + \j/ sin<9 w + (<p + \jf cos 6) z 
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V- Référentiel principal d’inertie 


V-l. Constatation 

Un référentiel est principal d’inertie s’il vérifie l’une des options suivantes : 
Option 1 : il contient l’axe de révolution du système. 

Option 2 : il contient deux axes de symétrie pour le système. 

Option 3 : il contient deux plans de symétrie pour le système. 


V-2. Conséquences 

Si St est un référentiel principal d’inertie alors la matrice d’inertie du système sera diagonale 


par rapport à St. 


C’est-à-dire : 


n G (s) 


A 0 0 
0 6 0 
0 0 C 


Remarques : 

> Si le système admet une symétrie sphérique alors A = B = C. 

y Si le système admet une symétrie cylindrique par rapport à Z (par exemple) alors A 

= B 

y n G (S) est la matrice centrale principale. 

> rus) = n 0 (S) + n. 4 (G,m) 

Pour déterminer U a (G ,m) , il suffit de déterminer le vecteur AG et l’exprimer dans la 
base où n G <s) est diagonale. Par exemple, si on trouve que : 


Alors : 


AG 


au + b w + c z 


(c’est le cas général) 


H A (G,m ) 


m(b 2 + c 2 ) 
-ni a b 
-ni a c 


— m a b 
m(a 2 + c 2 ) 
-m b c 


-m a c 
-m b c 


m(a 2 +b 2 ) 


(u,W,Z) 


Remarque : 


Si AG = c Z alors il suffit de donner à a et b la valeur 0, soit alors : 
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IIa(G,»î) 


me 2 0 0 

0 me 2 0 

0 0 0 --- 

- ( U,w,z ) 


VI- Energie cinétique 

• Si le système admet un point fixe A, par exemple, alors : 

T{S/R) = j Q(S/R 0 ) • n A (S) ■ Q(S/R 0 ) 

• Si le système n’admet pas de point fixe, alors : 

T(S/R 0 ) = ^mV 2 (G/R 0 ) + j Q(S / R 0 ) ■ U A (S) • Q(S / R 0 ) 


Remarque : 

Si le système Z est formé par (5^ ) U ( S 2 ) , alors : 

T(L/R 0 ) = T(S l /R 0 ) + T(S 2 /R 0 ) 

VII- Energie potentielle 

VII- 1. Energie de la pesanteur : 

Par définition U p (pesanteur) est définie par : 


U p (pesanteur) = - P ■ OG + este 



Centre du Centre d’inertie du 

référentiel fixe système 


Si l’axe Zq est orienté vers le bas alors P = m gz 0 , donc : 


U p (pesanteur) = -m g z 0 • OG + este 


Si l’axe Zq est orienté vers le haut alors P = -m gz 0 , donc : 


U p (pesanteur) = m g z 0 ■ OG + este 
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VII-2. Energie élastique : (énergie potentielle du ressort) 

Par définition U p (ressort) est : 


U (r es sort) = 


2 


K AX 2 


avec AX Ljj na j e L initiale 


Remarque : 

Si le système Z est formé par (5, ) U (S 2 ) , alors : 

u p (L) = u p (S l ) + u p (S 2 ) 

VIII- Equations de Lagrange 
VIII- 1. Système conservatif 

On dit qu’un système est conservatif si et seulement si les liaisons sont parfaites, 
indépendantes du temps et les forces qui travaillent dérivent d’une énergie potentielle. 

Remarque : 

On dit qu’une liaison est parfaite si et seulement si la puissance virtuelle est nulle. Cette 
puissance sera nulle si et seulement si : 

- Le mouvement s’effectue sans frottement, 

- Le point de la liaison est fixe, 

- Le roulement s’effectue sans glissement. 

VIII-2. Champ des déplacements virtuels 

On définit le déplacement virtuel d’un système comme étant un changement dans la 
configuration de ce système résultant d’un changement arbitraire très petit des coordonnées, 

noté S r L . compatible avec les forces et les liaisons imposées au système à l’instant donné t. 
Donc, ce déplacement virtuel n’implique que les déplacements des coordonnées et n’introduit 
aucune variation du temps St . 
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VIII-3. Equations de Lagrange 

Afin de former les équations de Lagrange, il faut établir une expression portant sur les 
déplacements virtuels des coordonnées généralisées qui sont indépendantes les unes des 
autres. 

Le passage des aux q, se fait à l’aide des équations de transformation : 

r i = r i (qvq 2 ’ -,q n ^) 


T} d r i 

Et aussi en utilisant les règles usuelles aux dérivées partielles. Ainsi la vitesse v ; = peut 
s’écrire à l’aide des q ,■ de la façon suivante : 


^ d n . d n 


De même, le déplacement virtuel arbitraire peut être associé aux déplacements virtuels par 


s “j 

j 


Le travail virtuel des forces Fi peut s’écrire, à l’aide des coordonnées généralisées, sous la 
forme suivante : 



J^FiSn = 


i i,j 

Avec : 

Qj , P- 

dqj 

Or: 

Fi = ntf Fj 



Xfiy 5( lj 


45 


i,j 1 7 


•J 


3(jf, 


d 2. dri 
2,-r t (m i — Sqj 
Tjdt dqj 


Or: 


3r; 


3 r i dt 


d 2 ri 
dt 


dF ; 


d V: 


dqj dt dq, d : q, dqj d<7 y 

a? 


Donc : 


dr 

ij "dj 


d 


dVi 


T,Fijr-*qj = Z*o»i^ a . 

dç, fj dt dq j 


— nt , V; 


dV l 
dq t 


)]<??/ 


= Six(^-4”< C » -4-4 m < r 2 )] s ii 

tj dt dqj 2 dqj 2 


Or: 


T- =-//i, V/ 


énergie cinétique de la particule i. 


^ ri ^ ^ 3r f 3r, 

i.j d( 'i U dt d( Ij d( ii 

^ r d ,dT^ d T is 

= 2-h-(vr-)-â — ]S Vj 

7 (lt dqj dqj 


On a: T F 'J^ S 9j = I^Qjdqj 

i-j d( lj j 

=> Z(4-(f^)-| = ° 

7 dt dqj dqj 

Comme est arbitraire alors : 

= o 

dt dqj dqj 1 

Dans le cas d’un système conservatif, les forces qui travaillent dérivent d’une énergie 
potentielle Up qui ne dépend pas du temps, et on peut écrire : 
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II 

<1 

U p donc Qj = 

= 

i 


D’où 

d ( dT ^ d T 
dt d qj d qj 

II 

1 

QJ OJ 

1 -ss C3 

ta 


=> 

d d 

(T U p )~ 

dt d qj 


= 0 


On pose L = T -U p le Lagrangien du système, on aura : 

= o 

dt d q j d qj 

VIII-4. Intégrale première de l’énergie 

Si le système est conservatif, alors on peut écrire l’intégrale première de l’énergie : 

T(S/R 0 ) + U p (S) = este 


VIII-5. Le Lagrangien du système 

Par définition : 


L(S/R 0 ) = T(S/R 0 ) - U p (S) 


VIII-6. Intégrale première du mouvement 
Définition : 

C’est le cas où on trouve qu’une fonction des paramètres du problème est égale à une 
constante 


IX- Positions d’équilibre 

Si le système est conservatif, on détermine les positions d’équilibre à l’aide de la résolution 

dU p 

du système d’équations — — = 0 où <?,- sont les paramètres principaux du système. 

d q t 


X- Etudes de la stabilité d’une position d’équilibre 

Si le système est conservatif, on peut déterminer si une position d’équilibre est stable à l’aide 
du théorème de Lejeune Dirichlet par : 


-On définit le Jacobien, 
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-On détermine la trace et le déterminant du Jacobien. 


Remarques : 

• Une position d’équilibre sera stable si, et seulement si, trace(J) > 0 et det (7) > 0 . 

• Si au moins l’un de ces deux quantités est <0 on ne peut pas conclure à l’aide de ce 
théorème. 


Exemple : 

Supposons que U p (S) = mglc osé? + mgRc os(p + este et que le système 
est conservatif. 

Position d’équilibre : 


su, 

de 


= — m g l sin 0 = 0 


0 = 0 ,71 


W P 

dtp 


= - m g R sin ç = 0 


(p = 0 , n 


Alors les positions d’équilibre sont : 


Ei (0i =0,çi = 0); E 2 (02 = 0,(p 2 =n); E 3 (0 3 = n, <pi= 0) ; E 4 (04 = n, ç 4 = n) 


Etude de la stabilité : 




= — m g l cos0 


= - m g R cosç> 



dOdtp 



dtpdO 


D’où : 


J = 


-mg IcosO 
0 


0 

-mg Rcostp 
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9 9 

Et: trace{J) = -m g l cosû -m g Rcosç et det(/) = m g IR cosOcosç) 

Pour la position 1 : 

trace(J ) = -m g l - mg R <0 

Donc : on ne peut pas conclure 
Pour la position 2 : 

detC J) = -m 2 g 2 IR < 0 

Donc : on ne peut pas conclure 
Pour la position 3 : 

detC J) = -m 2 g 2 IR <0 

Donc : on ne peut pas conclure 
Pour la position 4 : 

trace(J ) = m g l + m g R >0 
det(/)= m 1 g 1 IR > 0 

Donc : E4 est stable 
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